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ABSTRACT 



Aims. We study the radio emission mechanism of the FR-I AGN NGC 5128 (Centaurus A) 

Methods. We determine the centimeter and millimeter continuum spectrum of the whole Centaurus A radio source and 

measure at frequencies between 86 GHz (3.5 mm) and 345 GHz (0.85 mm) the continuum emission from the active 

radio galaxy nucleus at various times between 1989 and 2005. 

Results. The integral radio source spectrum becomes steeper at frequencies above 5 GHz, where the spectral index 

changes from aiow = -0.70 to Uhigh = -0.82. The SW outer lobe has a steeper spectrum than the NE middle and outer 

lobes {a = -1.0 vs -0.6). Millimeter emission from the core of Centaurus A is variable, a variability that correlates 

appreciably better with the 20-200 keV X-ray variability than with 2-10 keV variability. 

Conclusions. In its quiescent state, the core has a spectral index a — -0.3, which steepens when the core brightens. The 

variability appears to be mostly associated with the inner nuclear jet components that have been detected in VLBI 

measurements. The densest nuclear components are optically thick below 45-80 GHz. 

Key words. Galaxies - individual: NGC 5128 - centers; radio continuum: Centaurus A - AGN 



1. Introduction 

Centaurus A (NGC 5128) is the nearest (D = 3.4 Mpc) 
radio galaxy with an active nucleus. Because of this, it 
has been extensively studied at various wavelengths over 
a very wide range in the electromagnetic spectrum; much 
of the observational database was reviewed by Israel (1998). 
The radio source spans more than eight degrees on the 
sky and shows structure on all scales down to milliarc- 
seconds where VLBI techniques (e.g. Horiuchi et al. 2006) 
reveal an extremely compact core and bright nuclear jets. 
Notwithstanding the prominence of Centaurus A, relatively 
few measurements of the total radio emission of the source 
exist. The first comprehensive study of its decimetric emis- 
sion and polarization by Cooper et al (1965) has, in more 
recent times, been supplemented only by the measurements 
summarized by Alvarez et al. In particular, no measure- 
ments of the Centaurus A total flux density at frequencies 
above 5 GHz are found in the astronomical literature. In 
this paper, we extend that range by providing new deter- 
minations between 20 GHz and 100 GHz derived from the 
WMAP all-sky surveys. 

On much smaller scales, the flux densities of the com- 
pact Centaurus A core are also of interest. The central re- 
gion of Centaurus A has been imaged at very high res- 
olution with (space) VLBI techniques (Meier et al. 1989; 
Tingay et al. 1998, Fujisawa et al, 2000, Tingay, Preston 
& Jauncey 2001, Horiuchi et al. 2006) up to frequencies 
of 22 GHz, where both the selfabsorbed synchrotron emis- 
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sion from the core and the free-free absorbing foreground 
start to become optically thin (Tingay & Murphy 2001). 
Although Kellerman, Zensus & Cohen (1997) used the 
VLBA to observe Centaurus A at the signiflcantly higher 
frequency 43 GHz, they did not actually measure nuclear 
flux densities. Hawarden et al. (1993) have provided the 
most recent systematic attempt to derive the core spectrum 
from single-dish measurements, and they conclude that it 
is essentially flat at millimeter wavelengths, a conclusion 
reiterated by Leeuw et al (2002). 

Variability of the emission at frequencies of 22 and 43 
GHz was explicitly noticed by Botti & Abraham (1993), 
de Mello & Abraham (1990), and Abraham et al. (2007). 
However, these results, as well as those by Tateyama & 
Strauss (1992), generally suffered from low spatial reso- 
lutions of 2-4 arcmin that do not allow separate identi- 
fication of emission from the nucleus, the jets, and their 
surroundings. In this paper, we provide new flux density 
measurements of the compact core in the 90-345 GHz 
(sub) millimeter range, at resolutions between 0.8 and 0.2 
arcmin. As these measurements were collected over a time 
period spanning more than a decade, we are also in a posi- 
tion to study the variable behavior of this emission. 

2. Observations and reduction 

2.1. WMAP data 

The WMAP mission and its data products have been de- 
scribed in detail by Bennett et al (2003a, b, c). In our anal- 
ysis, we used the data products from the official WMAP 
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Fig. 1. Maps of the radio continuum emission of Centaurus A at (left to right) 23, 33, 41 GHz, 61 GHz, and 93 GHz. 
All images are at the nominal WMAP resolution. Contour levels are drawn at (23 GHz) 0.1, 0.19, 0.38, 0.74, 1.4, 2.8, 
5.5, 10.7 mK, (33 GHz) 0.1, 0.19, 0.36, 0.68, 1.3, 2.4, 4.6 mK; (41 GHz) 0.1, 0.19, 0.34, 0.64, 1.2, 2.2, 4.1, 7.6 mK; (61 
GHz) 0.1, 0.17, 0.31, 0.53, 0.94, 1.6, 2.9, 5.0 mK; (93 GHz) 0.1, 0.16, 0.25, 0.40 0.64 1.0, 1.6, 2.6 mK. Note the persistent 
diffuse emission features at +1.5°, -4° and at +2°, +8°, which are probably unrelated Galactic foreground emission (see 
Sect. 2.1.). 
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Fig. 2. Polarized vectors of the Centaurus A emission at (left to right) 23, 33, 41 GHz, and 61 GHz; the length of the 
vectors is proportional to the polarized intensity, but scaling is arbitrary. Far right: masks used to determine flux densities 
in Table[TJ 

Table 1. WMAP observations of Centaurus A 

Frequency Flux Densities Polarized Flux Densities 

Total NE Lobes+Core SW Lobes Total NE Lobes+Core SW Lobes 
(GHz) Jy 



23 


112±13 


76±9 


29±4 


6.8±0.5 


5.0±0.3 


1.8±0.2 


33 


82±12 


62±8 


20±4 


4.6±0.6 


3.6±0.3 


0.8±0.2 


41 


72±9 


56±6 


16±3 


3.8±0.5 


3.4±0.5 


0.4±0.2 


61 


56±15 


48±12 


8±3 


4.0±1 


4.0±1 


— 


93 


41±18 


47±18 


-6±9 


3.5±1 


3.5±1 


— 



5-year release, specifically the fully reduced and calibrated 
Stokes I, Q, and U maps of the entire sky, in nested 
HEALPIX format (Hinshaw et al. 2008). The maps were 
observed at frequencies (assuming synchrotron spectra) v 
= 22.5, 32.7, 40.6, 60.7, and 93.1 GHz with resolutions of 
53, 40, 31, 21, and 13 arcmin respectively. As a first step, 
the HEALPIX data maps were converted with the use of 
standard software to flat maps in Zenithal Equal Area pro- 
jection with pixel solid angles of p(o) = 1.90644 x 10^^ 
and intensities in mK. We first attempted to determine inte- 
grated flux densities by fitting in each map two-dimensional 
Gaussians to the emission corresponding to Centaurus A. 
However, with increasing frequency, both the beamsize and 
the total flux density decrease, resulting in a steep decline in 
flux density per beam. Our fitting results, not surprisingly, 
suggests that we were systematically missing low-surface- 
brightness extended emission. For this reason, we also inte- 
grated over the full extent of the Centaurus A radio source 
as shown in the maps by Cooper et al. (1965) - the mask 



used is depicted in Fig. [21 bottom right. We derived the flux 
densities given in Table[T]from the summed values with the 
conversion factor Jy/mK = 30.7 x p(o) x v^ . This yielded 
the flux densities given in Tablc[l] which correspond to the 
maps shown in Fig. [21 We verified the accuracy of this pro- 
cedure by also applying it to the strongest WMAP point 
sources J0319-I-4131, J0322-3711, J0423-0120, J1229-0203, 
J1230+1223, and J1256-0547 (Wright et al. 2008). For po- 
larization measurements, the Q and U maps were reduced 
using the same conversion factors as the I maps; the po- 
larized intensity was computed as P = a/Q^ + C/^. Only 
at 23 GHz are signal-to-noise ratios high enough to yield 
practically identical results for the gaussian-fitting and inte- 
gration methods. We consider the integrated flux densities 
at 23, 33, and 41 GHz to be quite reliable, whereas at 61 
GHz and 93 GHz much of the integration was over areas 
with a surface brightness very close to the noise, especially 
in the southern lobe. 
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In the total-intensity (I) maps, the spatial resolution 
blends all components into the dominant emission of the 
NE and SW radio lobes, respectively. Even at the highest 
frequency, the compact core cannot be separated from the 
intense NE lobe. This is different in the polarized intensity 
(P) maps. In the 23, 33, and 41 GHz P-maps, the com- 
pact core can be seen separately. At the higher frequencies, 
the signal-to-noise ratio of the P-maps is too low to yield 
useful results. We note that the central 22 GHz flux den- 
sities of Cen A typically range between 16 and 29 Jy in a 
4' beam (Fogarty & Schuch 1975; Botti & Abraham 1993), 
i.e. roughly 15-25% of the 23 GHz WMAP flux densities. 
At the same frequency, the VLBI milli-arcsec core (Tingay 
et al. 2001) contributes only 2% to the measured WMAP 
emission. 

Figure[2] shows the presence of an extended diffuse emis- 
sion feature to the east of the SW lobe. This feature can be 
seen at all frequencies and is undoubtedly real. However, 
we do not consider it to be related to Centaurus A. (1) 
Large-scale maps show significant foreground emission from 
the Milky Way reaching up to the Galactic latitude of 
Centaurus A (see Bennett et al. 2003b). (2) Inclusion of 
this feature in the flux densities of the SW lobe causes its 
emission to first decrease with frequency but then to make 
a surprising jump at 61 GHz, and at 93 GHz the weak 
SW lobe would have become as strong as the NE lobe -f 
core emission. (3) The lobes of Centaurus A are clearly 
polarized, whereas the diffuse feature does not show polar- 
ization at any frequency (Fig.[2l bottom row), which also 
argues against the feature being an extension of the SW 
lobe. Finally, (4) there is a similar but more obviously un- 
related feature at +2°, 8°. 

2.2. SEST data 

The SEST was a 15 m radio telescope located at La Silla, 
Chile, [j The main purpose of our observing program with 
the SEST was to obtain high-resolution molecular-line ab- 
sorption spectra against the continuum of the Centaurus 
A core, with integration times ranging from 30 minutes to 
4 hours per molecular transition per session, obtained by 
co-adding individual samples usually of 5 minutes each. 

In the time period 1989-1993, observations were made 
with Schottky barrier diode receivers providing single- 
sideband (SSB) system temperatures, including the sky, of 
500 K (3 mm band) and 1000 K (1.3 mm band). From 1996 
onwards until the closure of the SEST in 2003, observations 
were made with Superconductor-Insulator-Superconductor 
(SIS) tunnel junction receivers operating simultaneously in 
either (i) the 3 mm and 2 mm bands, or (ii) the 3 mm and 
1.3 mm bands. Typical system temperatures, including the 
sky, were 220 K in the 3 mm band and 330 K in the 1.3 mm 
bands. Continuum pointings were frequently performed on 
Centaurus A itself; monitoring the continuum level during 
the observations provided an additional and very valuable 
check on tracking accuracy. During the long integrations, 
the r.m.s. pointing was on average on the order of 5" or bet- 
ter. All observations were made in double-beamswitching 
mode, with a 6 Hz switching frequency and a throw of 12' 
because this procedure yielded excellent and stable base- 



lines. SEST beam sizes and main-beam efficiencies at the 
various observing frequencies were listed by Israel (1992); 
they were 41"-57" and 0.70-0.75 in the 3 mm band, and 23" 
and 0.5 in the 1.3 mm band. Whenever feasible, we used 
the observatory high- and low- resolution acousto-optical 
spectrometers (AGS) in parallel. 

For this paper, we extracted the continuum flux den- 
sity from measurements with the low-resolution backcnd 
(bandwidth 500 MHz, resolution 1 MHz). To this end, we 
averaged, for each 5-minute sample, the continuum levels in 
the 0-300 km s""'^ and 800-1100 km s^^ velocity intervals, 
straddling but avoiding the central line emission proflles. 
Individual sample continuum antenna temperatures were 
then converted to flux densities. The conversion typically 
ranged from 25 Jy/K at 86 GHz, 27 Jy/K at 115 GHz, 
30 Jy/K at 147 GHz, to 41 Jy/K at 230 GHz. The flux 
densities listed in Tabled for each observed frequency are 
the means of the individual samples observed at the fre- 
quency and run listed. Observing runs typically lasted four 
to eight days; the median date is given. In determining the 
means, we left out samples negatively affected by pointing 
drift, which were easily recognized by a slow but steadily 
decreasing continuum level. 

From 1996 to 1999, the Centaurus A continuum inten- 
sity was also monitored every few months at frequencies 
of 90, 150, and 230 GHz in a dedicated program. The ob- 
servations were performed in the dual beam switch mode 
with a 12' throw, and the data were taken simultaneously 
at either 90/150 or 90/230 GHz. Individual exposures were 
kept short to minimize the effect of the atmosphere and 
to allow a flltering of spectra with poor quality. The r.m.s. 
values from the pointing runs were typically better than 3". 
For each of the individual spectra, the continuum flux level 
was determined by fltting a zeroth-order line. Finally, the 
average was determined by a least-square fit to the fitted 
continuum levels and error of all the individual measure- 
ments. 

The quality of the data is a function of the stability 
of the sky (which determines the accuracy of the double- 
beamswitch signal subtractions) and the accuracy of the 
pointing. As the pointing depends on the proper evaluation 
of a series of boresight observations taken in sequence, it is 
itself a function of sky stability. In general, observations in 
the 3 and 2 mm windows profited from transparent sky con- 
ditions superior to those encountered in the 1 mm window. 
In addition, a particular r.m.s. pointing error has a greater 
effect on smaller beamsizes, i.e. at higher frequencies. For 
this reason, observations around e.g. 90 GHz have the high- 
est quality, and those at 220/230 GHz are associated with 
larger relative errors. 

2.3. JCMT data 

From 2002 through 2005, similar molecular line obser- 
vations were made with the 15m James Clerk Maxwell 
Telescope (JCMT) on Mauna Kea (Hawaii). 0. They were 
made in beamswitching mode with a throw of 3' in azimuth 
using the DAS digital autocorrelator system as a backend. 



The Swedish-ESO Submillimeter Telescope was operated 
jointly by the European Southern Observatory (ESO) and the 
Swedish Science Research Council (NFR). 



The James Clerk Maxwell Telescope is operated on a joint 
basis by the United Kingdom Particle Physics and Astrophysics 
Council (PPARC), the Netherlands Organisation for Scientific 
Research (NWO), and the National Research Council of Canada 

(NRC). 
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Table 2. SEST observations of the Ccntaurus A nucleus 



Day 


Date 


Frequency 


Brightness 


Flux Density 


Day 


Date 


Frequency 


Brightness 


Flux Density 






V (GHz) 


Tn,b (mK) 


S. (Jy) 








V (GHz) 


T^b (mK) 


S. (Jy) 


1989.29 


Apr 15 


87 


373±45 


7.2±0.9 


1998.61 


Aug 


15 


90 


286±18 


7.1±0.7 


1989.50 


Jul 3 


89 


391±15 


7.3±0.7 








230 


138±10 


5.7±0.6 






98 


393±7 


7.5±0.8 


1998.99 


Dec 


30 


90 


321±13 


8.0±0.8 






110 


313±20 


5.9±0.6 








150 


226±16 


6.8±0.7 






113 


353±9 


6.6±0.7 


1999.10 


Feb 04 


90 


277±16 


6.9±0.7 






220 


223±24 


5.0±0.7 








150 


184±20 


5.5±0.6 






227 


208±44 


4.3±0.8 


1999.12 


Feb 11 


90 


304±16 


7.6±0.8 


1991.03 


Nov 01 


88 


371±20 


7.1±0.8 








150 


283±37 


8.5±1.2 






89 


369±21 


7.0±0.8 


1999.13 


Feb 18 


90 


333±24 


8.3±0.6 






91 


361±20 


6.9±0.8 








150 


282±33 


8.5±1.0 






114 


361±25 


6.8±0.9 


1999.38 


May 


17 


86 


372±45 


7.0±0.9 


1992.53 


Jul 13 


89 


516±11 


9.7±1.0 








90 


373±45 


7.1±0.9 






91 


516±7 


9.7±1.0 








98 


362±45 


6.9±0.7 






99 


488±31 


9.2±1.2 








110 


377±45 


7.2±0.7 






218 


263±11 


5.5±0.6 








169 


168±45 


5.5±0.6 






220 


326±60 


6.8±1.3 








226 


276±45 


5.5±0.6 






226 


289±35 


5.8±0.8 








230 


250±45 


5.1±0.5 






268 


210±13 


4.8±0.7 








267 


190±45 


4.4±0.8 


1993.55 


Jul 20 


218 


384±40 


8.0±0.8 


1999.49 


Jun 


26 


90 


394±14 


9.8±1.0 






226 


354±40 


7.3±0.7 








150 


297±16 


8.9±0.9 






268 


329±37 


7.6±0.8 


1999.61 


Aug 


12 


90 


299±18 


7.5±0.8 






346 


292±30 


7.6±0.8 








115 


285±40 


7.7±1.1 






357 


238±35 


5.6±0.6 








150 


299±18 


9.0±0.9 


1996.67 


Aug 31 


110 


353±17 


6.7±0.7 


1999.98 


Dec 


25 


90 


292±25 


7.3±0.7 






115 


333±31 


6.3±0.6 








150 


250±26 


7.5±0.8 






145 


275±14 


5.5±0.6 


2000.37 


May 


15 


85 


408±21 


7.8±0.8 


1996.97 


Dec 19 


90 


410±20 


10.3±1.0 








89 


397±19 


7.4±0.8 






229 


200±20 


8.2±0.8 








98 


374±24 


7.2±0.7 


1996.99 


Dec 29 


90 


380±20 


9.5±1.0 








111 


374±22 


7.1±0.7 






229 


160±20 


6.6±1.6 








219 


233±20 


4.8±0.5 


1997.10 


Feb 06 


90 


350±10 


8.8±0.9 








262 


205±68 


4.7±0.5 






150 


250±20 


7.5±0.8 








267 


227±30 


5.2±0.6 


1997.15 


Feb 21 


90 


370±10 


9.3±0.9 


2002.07 


Jan 


25 


87 


521±30 


9.8±1.0 






150 


290±10 


8.7±0.9 








91 


495±25 


9.4±0.9 


1997.30 


Apr 19 


90 


360±30 


9.0±0.9 








109 


406±50 


8.6±1.0 


1997.68 


Sep 05 


87 


440±20 


8.3±0.8 








267 


243±70 


5.5±1.6 






89 


420±50 


7.9±0.9 


2002.36 


May 


10 


84 


482±11 


8.9±0.9 






113 


353±55 


6.6±0.8 








85 


477±5 


9.0±0.9 






145 


318±30 


6.3±0.6 








87 


475±27 


8.9±0.9 






150 


393±30 


6.0±0.6 








88 


459±23 


8.6±0.9 






220 


265±50 


5.4±1.0 








89 


476±15 


8.9±0.9 


1998.49 


Jun 28 


90 


260±20 


6.5±0.7 








93 


511±14 


9.5±1.0 






150 


300±20 


9.0±0.9 








97 


461±15 


8.8±0.9 


1998.50 


Jul 03 


89 


424±20 


8.0±0.8 








109 


438±25 


8,3±0.8 






98 


370±20 


7.0±0.7 








147 


375±17 


7.4±0.7 






218 


335±20 


6.9±0.7 








218 


213±17 


4.4±0.5 






227 


212±20 


4.0±0.8 








220 


300±33 


6.1±0.6 


1998.59 


Aug 08 


90 


290±10 


7.3±0.7 








227 


290±14 


5.9±0.6 






150 


240±10 


7.2±0.7 








230 


280±28 


5.7±0.6 


1998.60 


Aug 14 


90 


280±20 


7.0±0.7 








262 


204±42 


4.6±0.7 






230 


140±20 


5.7±0.9 


2002.96 


Apr 15 


89 


384±43 


7.2±0.8 


















220 


227±45 


4.6±0.9 


2003.16 


Feb 26 


85 


421±22 


8.0±0.8 


2003.21 


Mar 16 


83 


478±25 


9.1±0.9 






87 


437±13 


8.2±0.8 








85 


455±16 


8.7±0.9 






88 


420±16 


7.9±0.8 








86 


532±20 


lO.Oil.O 






89 


416±24 


7.8±0.8 








89 


508±41 


9.5±1.0 






107 


432±17 


8.1±0.8 








91 


479±29 


8.9±0.9 






218 


290±23 


6.0±0.6 








104 


467±17 


8.8±0.9 






220 


278±45 


5.7±0.9 








110 


438±34 


8.3±0.8 






227 


286±32 


5.8±0.7 








169 


353±19 


7.2±0.7 






262 


241±38 


5.4±0.8 








220 

226 

227 


220±30 
298±28 
244±30 


5.6±0.9 
6.0±0.6 
6.1±0.7 
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Table 3. JCMT observations of the Centaurus A nucleus 



Day Date Frequency Main-Beam Flux Density Day 

V Brightness 8,^ 

(GHz) T^b (mK) (Jy) 



Date Frequency Main-Beam Flux Density 
V Brightness Si/ 

(GHz) T^b (mK) (Jy) 



1997.09 


Feb 01 


330 


2002.28 


Apr 12 


230 
268 


2002.47 




268 


2003.24 


Mar 30 


268 
330 
345 



327±22 



5.8±0.6 


2003.36 


May 05 


265 


6.3±0.9 


2003.50 


Jul 01 


265 


5.9±0.9 






330 


6.2±0.9 


2004.30 


Apr 17 


345 


6.1±0.9 


2004.57 


Jul 26 


330 


7.9±1.2 


2005.23 


Mar 23 


330 


5.2±0.8 









5.4±0.7 
4.1±0.6 
8.4±1.3 
7.2±1.1 
6.3±0.9 
4.4±0.7 




0.1 1 

Freqency (GHz) 



Fig. 3. The area-integrated continuum spectrum of the 
Centaurus A radio source from meters to millimeters. Filled 
circles represent integrated flux densities from the literature 
(cf. Alvarez et al. 2000) and from Table[TJ For comparison, 
open triangles represent integrated flux densities including 
the SE feature which we suspect to be unrelated Milky Way 
foreground (see Sect. 2.1) whereas open circles represent the 
flux densities obtained from Gaussian fitting of the WMAP 
emission which is biased to regions of relatively high sur- 
face brightness. Solid lines mark fits to the literature and 
WMAP data. 



As the low declination of Centaurus A brings the source to 
less than 30° above the horizon as seen from Hawaii, ob- 
servations were performed in timeslots covering one hour 
at either side of transit. For all spectra, we scaled the ob- 
served continuum antenna temperature to flux densities by 
assuming aperture efficiencies ?7ap = 0.57 at 230-270 GHz 
and 77mb = 0.49 at 330-345 GHz. 

3. Results and analysis 

3.1. The overall Centaurus A radio source spectrum 

In Fig.[3]we show the spatially integrated radio spectrum of 
Centaurus A from 85 MHz to 90 GHz. The lower- frequency 
flux densities have all been taken from the literature as com- 
piled, corrected, and discussed by Alvarez et al. (2000), who 
found a spectral index a = -0.70±0.01 (5^ oc ly") corre- 



sponding to the solid line in Fig.[3l The higher-frequency 
points were taken from Table[T] A least-square fit indi- 
cates a slightly steeper spectrum with spectral index a 
= -0.82±0.07 and a spectral turnover by Aa w —0.12 
somewhere between 5 and 20 GHz (see also Junkes et al. 
1997). The more poorly defined spectrum of the SW lobe 
is steeper than that of the NE lobe with asw = — 1.0 vs 
CKNE-i-c = —0.6, but this difference is barely significant. 

In these spectra and flux densities, we do not include the 
contribution by the extended emission feature east of the 
SW lobe discussed above. If we were to do so, its relatively 
strong presence at the higher frequencies would flatten the 
high-frequency spectrum to the one indicated by the open 
triangles in Fig. [3] (with a = —0.64 ± 0.09). The slope 
of that spectrum is very similar to that at the lower fre- 
quencies, but the spectrum itself is offset in flux density by 
about a factor of two and does not fit very well onto the 
lower-frequency spectrum. It is clear from Fig. [3] that the 
emission from the diffuse feature becomes negligible com- 
pared to the Centaurus A total at frequencies below fre- 
quencies 20 GHz. The low Galactic latitude of Centaurus 
A also has some effect at much lower frequencies. For in- 
stance. Cooper et al. (1965) briefly discuss the presence of 
a relatively weak spur of Galactic foreground synchrotron 
emission at the southeastern tip of Centaurus A. The lo- 
cation and extent of this spur is indicated in their Fig. 1. 
The WMAP diffuse feature is close to, but not coincident 
with, the spur and appears to be unrelated to it. Indeed, 
none of the radio continuum maps shown by Cooper et al. 
(1965) and Junkes et al. (1993) in the range of 0.4 to 4.8 
GHz shows emission coincident with the WMAP feature. 

There is some evidence that the ratio i?NS of the flux 
densities of the NE lobes (including the central core emis- 
sion) and the SW lobes increases with frequency. At 0.4 
and 1.4 GHz, i?NS = 1-5 (Cooper et al. 1965, Junkes et 
al. 1993), but this has increased to i?NS ~ 3 in the 23-41 
GHz range (Table[T]). Even if we were to add the suspected 
foreground emission to that of the SW lobe, we still would 
have a ratio i?NS = 2 in the well-determined 23-43 GHz 
range. We note that Alvarez et al. (2000) have found that, 
over the full 0.08-43 GHz range, at least the inner lobes 
also have a ratio i?NS = 1.5. For the giant owier lobes they 
find a lower ratio _Rns = 0.9, but this only applies to the 
limited 0.4-5.0 GHz range. The meaning of these apparent 
changes in ratio is not clear to us. 

The fraction of polarized emission is well-defined be- 
tween 23 GHz and 41 GHz. The lower-frequency, higher- 
resolution maps by Cooper et al. (1965) and Junkes et al. 
(1993) show overall polarization increasing from P/I « 
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Table 4. 

cleus 



Normalized flux densities of the Centaurus A nu- 




Fig. 4. Flux densities of the Centaurus A nucleus at 90 
GHz in 1988-2005 taken from this paper (filled circles) and 
from published literature (open circles; see text). Dashed 
lines represent X-ray fluxes monitored by the BATSE and 
RXTE experiments on an arbitrary flux density scale 



12% at 960 MHz to P/I « 30% at 5 GHz/. We flnd much 
lower polarized fractions. The polarization of the entire 
source is within the errors constant over the 23-61 GHz 
range with P/I = 6.0 ± 0.3%. As Fig. [2] shows, the beam- 
averaged polarization direction is different for all three com- 
ponents. The low polarization fraction found by us most 
likely reflects the significant beam depolarization caused 
by the low resolution of the WMAP data. 

3.2. Millimeter-wave emission of the Centaurus A core region 

Although the extended radio emission from Centaurus A 
even at millimeter wavelengths significantly exceeds that 
of the core region, its surface brightness is low and de- 
clines rapidly with increasing frequency. Can it, however, be 
safely ignored in the flux density measurements of the core 
with the SEST beam? The highest frequency measured with 
WMAP, 93 GHz, corresponding to a wavelength of 3.2 mm, 
is in the lowest band that the SEST offered for ground based 
millimeter- wave molecular line observations. The measured 
gaussian-fitted fiux density at this frequency is S93(tot) = 
24±5 Jy, well above the flux density S = 9±1 Jy that single- 
dish measurements suggest for the Centaurus A core region. 
Excluding the core, the mean surface brightness at 90 GHz 
of the extended synchrotron emission component is about 
115 mJy/arcmin^. The contribution of the extended non- 
thermal component can thus safely be neglected in the rel- 
atively high-resolution observations of the core. 

At very high frequencies, thermal emission from dust 
may also contribute to measured intensities, but this is not 
the case in the frequency range of interest to us. The con- 
tinuum emission seen at resolutions of a few tenths of an 
arcminute by telescopes such as the SEST and the JCMT, 
effectively originates in a point source as can clearly be 
seen in the 850/im (345 GHz) SCUBA images published by 
Mirabel et al. (1999) and Leeuw et al. (2002). From these 
and our own data, we deduce that the extended emission, 
after correction for 20% CO line contamination, contributes 
less than 4% to the central peak intensity at 345 GHz, and 
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Index 
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— 





— 
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1989.50 
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1991.03 
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1992.53 
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5.7 


— 


-0.58 


1993.55 
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6.6 
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1996.67 
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5.4 
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-0.32 


1997.15 


9.3 


8.7 
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7.1 
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5.7 
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— 


1999.38 


7.2 


6.1 


5.3 
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— 
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— 


— 


2004.30 


— 


— 


— 


7.2 
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— 


— 


— 


4.3 


— 


Means 


8.1 


7.1 


5.9 


(6.1) 


-0.33 



this contribution drops rapidly with decreasing frequency. 
In the next section we will return to this issue. 

To study the temporal behavior of the mm-wave con- 
tinuum spectrum of the core, we have reduced the mea- 
surements listed in Table[3] to standard frequencies of 90, 
150, and 230 GHz. We performed a linear regression on all 
flux density/frequency pairs measured during a run, and for 
each run determined the best-fitting 'standard' fiux densi- 
ties by interpolation, as well as the effective spectral index 
between 90 and 230 GHz. As a consequence of this pro- 
cedure, the results for each run, given in Tabled are more 
robust than the individual measurements taken during that 
run. Typical uncertainties are about 0.3 Jy at 90 GHz and 
0.6 at 230 GHz. Although the sampling in time was not uni- 
form, we have given the mean fiux density over the whole 
period for each frequency at the bottom of Tabled) If more 
detailed information is absent, these values are the best 
guess for the core flux densities at any time not covered 
by the monitor. We have also listed the rather flat spectral 
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Fig. 5. Relative changes in 90 GHz flux densities as a func- 
tion of observing interval. Changes above the horizontal 
thin line are highly significant and not influenced by ob- 
serving errors or bias. 
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Fig. 6. Spectral index 90-230 GHz as a function of 90 GHz 
flux density. 



index defined by these 'best' points, with the important 
caveat that the three standard flux densities defining the 
means were not sampled at the same times. 



3.3. Core variability at mm wavelengttis 

In Fig.|3]we have plotted all observed 90 GHz flux densities 
from Table[2] and added the available (SEST) values from 
the literature (Tornikoski et al., 1996; Wiklind & Combes 
1997) fortunately filling in much of the gap created by our 
lack of monitoring Centaurus A in the period 1993-1997. 
It was clear from the start that the core of Centaurus A is 
variable. In Fig. [4] we illustrate this variability of core flux 
density. Over the period monitored, the flux density never 
drops much below 7 Jy, but at times may be higher than 
10 Jy. Sustained periods of relatively high flux densities ap- 
pear to have occurred in 1993/1994, 1997, and 2001/2002. 
There is no doubt that this variability is real. The ex- 
cursions exceed the uncertainty associated with individual 
plotted values, and the behavior is similar at all observed 
frequencies. At 90 GHz, the core flux density is seen to 
change by as much as 20% in a few weeks. This is illus- 
trated in Fig. [51 which plots the percentage change from 
one observing date to another as a function of the elapsed 
time. Over the observed period, the greatest changes (more 
than 25%) take place over time intervals typically between 
2 and 4 months. The high-frequency variability observed 
here should be compared to the variability deduced from 
the lower frequency monitoring observations by Botti & 
Abraham (1993) and Abraham (1996). Although there is, 
unfortunately, only limited temporal overlap, the pattern 
is very similar, with signiflcant and sometimes rapid inten- 
sity changes. In particular, enhanced radio emission is seen 
in 1992 in both our sample and in the data by Abraham 
(1996). The 43 GHz variability has a large amplitude of 8 - 
10 Jy (corresponding to factors of two change in brightness) 
and is about half as much (factor of 1.25) at 22 GHz. We 
note that those flux densities were determined with rela- 
tively large beams of 2.1' and 4.2', allowing signiflcant con- 
tributions to the emission especially at 22 GHz, not related 
to the core. 

The X-ray variability of the Centaurus A nucleus was 
monitored over the same period in the 20-200 keV band by 
the BATSE experiment onboard the Compton Gamma Ray 
Observatory from 1991 through 2000 and in the 2-10 keV 
band by the RXTE AU Sky Monitor from 1996 onwards. 
We took the public data from the NASA Goddard SFC and 
MIT websites and constructed averages over 15-day bins. 
The results are plotted with the 90 GHz flux densities in 
Fig.m The correspondence between the SEST 90 GHz and 
the BATSE X-ray flux densities is quite good, suggesting 
closely related, perhaps even identical, radiative processes 
. The correlation between the millimeter-wave continuum 
and the relatively soft X ray emission sampled by the RXTE 
experiment is much poorer. Minor X-ray excursions (on the 
order of a factor of two) occurred in 1991/1992, during 1994, 
in late 1996 and early 1997, during 1999, and 2001. They 
correlate with radio maxima, although unfortunately we are 
lacking millimeter continuum coverage for 2001. 

Although the behavior is the same at all frequencies ob- 
served by us, flux densities vary less at higher frequencies. 
As a result, the core spectral index changes with the degree 
of activity. In Tabled we have calculated spectral index val- 
ues a [Sy oc u") for the frequency range of 90 GHz to 230 
GHz. To guarantee a reasonable accuracy (on the order of 
0.10 or better in the index), we included only epochs at 
which measurements were available in three different fre- 
quency windows. In Fig. [5] we have plotted these spectral 
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indexes as a function of the 90 GHz flux density. It is quite 
obvious that the nim-wave spectrum of the Centaurus A 
core becomes steeper (smaller spectral index) when the ac- 
tual flux increases; i.e. increases are less in both absolute 
and relative senses at the higher frequencies. The points in 
Fig.[6l with their uncertainties of Aa w 0.10, can be fitted 
with a linear-regression line a = -0.087 Sgo + 0.311. For 
the observed lower limit of 7 Jy to the Centaurus A core 
flux density, this fit implies a fairly flat spectral index a = 
-0.3 (although the two actually measured values imply an 
even flatter spectrum with a — -0.22). Spectral indices a 
= -0.6 are approached at enhanced flux densities of 10 Jy. 
Variations in the spectral index between 22 and 43 GHz 
were likewise deduced by Botti & Abraham (1993) and plot- 
ted as a function of time in their Figure 2. If we plot their 
data again, now with the spectral index as a function of the 
relatively strongly varying 43 GHz flux density, we flnd a 
behavior opposite to that deduced by us for the 90 - 230 
GHz range. With relatively little scatter, 022-43 = + 0.11 
5*43 - 2.2; i.e. the spectrum flattens when S43 increases! 
We do not have a full explanation for this, but some re- 
marks are in order. First, we note that the 43 GHz flux 
density appears to vary even more than the 90 GHz flux 
density. Although there is very little temporal overlap, this 
is nevertheless consistent with our conclusion that the spec- 
trum steepens when flux densities increase. Second, the 22 
GHz emission would be expected to vary even more, but 
obviously fails to do so. This could be explained by high 
optical depths of the variable emission component at this 
frequency. 

3.4. Nature of the core emission 

In the quiet state, the core of Centaurus A has a high- 
frequency spectral index atigh ~ -0.3, but it must become 
optically thick around 80 GHz to allow flux densities to be a 
factor of two lower at 43 GHz than at 90 GHz. This implies 
that the emission at 22 GHz must arise almost entirely 
from non-core components. In the period covered by our 
observations, Centaurus A has not shown much activity, in 
any case less than in earlier decades (cf Abrahamfe Botti 
1993). Nevertheless, we did observe a few episodes with 
slightly enhanced activity. In these mildly 'active' states, 
flux densities increase, the high-frequency spectral index 
steepens to auigh — -0.6, and the spectral turnover shifts 
to lower frequencies of about 45-50 GHz. In this state, there 
is a measurable 2-4 Jy contribution from the core to the 
total 22 GHz emission from the central 4 arcmin, such as 
the 4 Jy seen in the VLBI measurements by Tingay et al. 
during the 1995 episode, but most of the 22 GHz emission 
still represents non-core components. 

Recent (April 2006) observations with the Smithsonian 
Milhmeter Array (SMA) at a resolution of 2" x 6" show a 
point source with a continuum flux density S23QGHZ — 6 
Jy (Espada et al. 2007). Although we do not have a single- 
dish measurement at the same epoch, a glance at Table|4] 
will show that this is, in effect, the total flux to be expected 
at this frequency. 

What does the observed core emission represent, and 
what is the meaning of the variability? Although our ob- 
serving beams were much larger, essentially all observed 
emission arises from a region with a size of a few arcsec at 
most, as indicated by the SMA observations just mentioned, 
as well as the measurements described below. VLBI images 



at frequencies of 5 and 8 GHz (Tingay et al. 1998, 2001; 
Tingay & Murphy 2001; Horiuchi et al. 2006) show very 
considerable source structure concentrated within about 60 
milliarcsec. We note that the submillimeter core itself suf- 
fers from both free-free absorption and synchrotron selfab- 
sorption at these frequencies (see e.g. Tingay & Murphy 
2001). Nevertheless, the images suggest that the emission 
observed by us arises from a nuclear source, as well as from 
bright nuclear jets extending over the better part of a par- 
sec. This complicates the interpretation of the observed 
variability, since we have no means of determining where 
precisely it originates. The time sequences registered by 
Tingay et al. (1998, 2001) strongly suggest that nuclear jet 
evolution is a major source of variability in the emission, 
as total 8 GHz flux densities vary between 4 and 10 Jy, 
Their maps also show movement and variability of indi- 
vidual components in the nuclear jets. The 8 GHz 'core' 
component varies much less (between 1.8 and 3.3 Jy with a 
median of 2.4 Jy over the period 1991-2000) and the 8 GHz 
and 22 GHz core flux densities are almost identical at the 
only two epochs (1995.88 and 1997.23) for which a reliable 
measurement exists (Tingay et al. 1998, 2001). However, in 
the much higher resolution (1.2x0.6 mas) maps by Horiuchi 
et al. (2006) this 'core' breaks up into at least four sepa- 
rate components. These extend over 7 milliarcsec (0.13 pc) 
and should presumably be considered as inner nuclear jet 
components. It is thus unclear whether the actual nuclear 
source (if any) itself is variable or not. 

Abraham et al. (2007) have interpreted rapid, day-to- 
day variations with amplitudes of 20% in the 43 GHz emis- 
sion observed over a three-month period in 2003 as evidence 
of free-free absorption of the nuclear source by clouds in 
the center of NGC 5128. It is, however, by no means clear 
that such an explanation could also apply to longer term 
variations, or indeed to the pattern of variability at higher 
frequencies. 

Figure[6] shows that the core spectrum steepens when 
it brightens, implying rapid energy loss of any electrons 
injected during active periods. We might assume that the 
emission from the Gen A core consists of a variable compo- 
nent and a component of constant emission. In that case, 
the spectral index of the variable emission is reasonably 
well-constrained with ayar = -1.6, especially if the 43 GHz 
peak excursions of about 9 Jy (Abraham 1996) correspond 
to the 3 Jy peak excursions at 90 GHz, which we do not 
know for certain. The relative flux densities of the two com- 
ponents and the spectral index of the constant component 
need an additional constraining assumption. The lowest 
flux densities measured at 43 GHz are 3-4 Jy, which puts an 
upper limit on the constant component. For instance, if we 
assume a flat spectrum for the latter {acst = 0.0), the mini- 
mum flux density of the variable component is S*^™ (var) = 
3 Jy. However, there does not seem to be a consistent model 
for the variation at all three frequencies without additional 
and ad- hoc assumptions. It appears that millimeter- wave 
VLBI monitoring is needed to provide the deflnitive an- 
swer to the question of the Centaurus A nuclear source 
properties. In particular, it would be quite interesting to 
compare the results of such monitoring with the results for 
the core of Centaurus A extracted by Meisenheimer et al. 
(2007) from mid-infrared interferometry at the ESO-VLT. 
After completion, the maximum resolution of the Atacama 
Large Millimeter Array (ALMA) will be about 50 mas at 
100 GHz. This will close the gap between the centimeter- 
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wavelength synthesis array resolutions and the centimeter- 
wavelength VLBI fields of view noted by Israel (1998, Sect. 
5.4), and thus be invaluable for studies of the evolution of 
the nuclear and inner jets. The resolution falls far short, 
however, of those needed to separate processes occurring in 
the nuclear jets and in the core components. If THz imag- 
ing is feasible at the longest baselines, the resulting ALMA 
resolution of 5 mas is better suited to such studies. ALMA 
sensitivities will be amply sufficient to show the core, but it 
remains to be seen how much of the subparsec jet structure 
will be discernible. 



4. Conclusions 

We have extended the spatially integrated continuum spec- 
trum of the Centaurus A radio source by using WMAP 
results between 23 and 93 GHz. The spectrum has now 
been explored over four frequency decades, from 10 MHz 
to 93 GHz. Between 10 MHz and 5 GHz, the spectral index 
a — —0.70, as determined in the literature. Longwards 
of 5 GHz the spectrum appears to steepen, with a = 
-0.83 ± 0.07. The flux density ratio Rm of the NE lobes 
(including the core) and the SW lobes increases with fre- 
quency from i?NS ~ 1-5 at 0.4-1.4 GHz to Rns ~ 3 in 
the 23-61 GHz range. 

We monitored the emission from the Centaurus A core 
component at frequencies between 90 and 230 GHz over 
more than ten years. The core is variable in emission, but 
over the period 1989-2003, the 90 GHz flux density from 
the compact core was never below 7 Jy. The spectrum of 
the quiescent core is relatively flat with a = -0.3. Variability 
was less at higher frequencies, implying steepening of the 
continuum spectrum simultaneous with core brightening. 
During active periods, the emitted flux increases, but the 
spectral turnover frequency decreases. The nuclear compo- 
nent is optically thick below 45-80 GHz. 

It appears that most if not all of the variability is as- 
sociated with the inner nuclear jet components that have 
been detected in VLBI measurements, but the mechanism 
of variability is not yet clear. 
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